Screening and identification of lncRNAs in preadipocyte differentiation in sheep.

Studies of preadipocyte differentiation and fat deposition in sheep have mainly focused on functional genes, and with no emphasis placed on the role that long non-coding RNAs (lncRNAs) may have on the activity of those genes. Here, the expression profile of lncRNAs in ovine preadipocyte differentiation was investigated and the differentially expressed lncRNAs were screened on day 0 (D0), day 2(D2) and day 8(D8) of ovine preadipocyte differentiation, with their target genes being predicted. The competing endogenous RNA (ceRNA) regulatory network was constructed by GO and KEGG enrichment analysis for functional annotation, and some differentially expressed lncRNAs were randomly selected to verify the RNA-Seq results by RT-qPCR. In the study, a total of 2517 novel lncRNAs and 3943 known lncRNAs were identified from ovine preadipocytes at the three stages of differentiation, with the highest proportion being intergenic lncRNAs. A total of 3455 lncRNAs were expressed at all three stages of preadipocyte differentiation, while 214, 226 and 228 lncRNAs were uniquely expressed at day 0, day 2 and day 8, respectively. By comparing the expression of the lncRNAs between the three stages of differentiation stages, a total of 405, 272 and 359 differentially expressed lncRNAs were found in D0-vs-D2, D0-vs-D8, and D2-vs-D8, respectively. Functional analysis revealed that the differentially expressed lncRNAs were enriched in signaling pathways related to ovine preadipocyte differentiation, such as mitogen-activated protein kinase (MAPK) pathway, the phosphoinositide 3-kinase protein kinase B (PI3K-Akt) pathway, and the transforming growth factor beta (TGF-ß) pathway. In summary, lncRNAs from preadipocytes at different stages of differentiation in sheep were identified and screened using RNA-Seq technology, and the regulatory mechanisms of lncRNAs in preadipocyte differentiation and lipid deposition were explored. This study provides a theoretical reference for revealing the roles of lncRNAs in ovine preadipocyte differentiation and also offers a theoretical basis for further understanding the regulatory mechanisms of ovine preadipocyte differentiation.

Spatiotemporal Expression and Haplotypes Identification of KRT84 Gene and Their Association with Wool Traits in Gansu Alpine Fine-Wool Sheep.

Keratin (K) is a major protein component of hair and is involved in hair growth and development. In this study, we analysed the expression, localization, and polymorphism of the K84 gene (KRT84) in Gansu Alpine Fine-wool sheep using immunofluorescence, RT-qPCR, and PARMS (penta-primer amplification refractory mutation system). Haplotypes of KRT84 were also constructed and their relationship with wool traits analysed. It was revealed that KRT84 was highly expressed in hair follicles, including the inner root sheath, outer root sheath, and hair medulla and at all six lamb ages investigated from 1 to 270 days of age. Three SNPs were detected in KRT84 exon 1, and they formed three haplotypes (named H1, H2, and H3) and six genotypes. Analyses revealed an association between haplotype combinations (diplotypes) and the mean fibre curvature, mean staple length, mean staple strength, mean fibre diameter, the coefficient of variation of fibre diameter, and comfort factor for these sheep. These results suggest that KRT84 is of importance in determining several key traits in Gansu Alpine Fine-wool sheep and that the gene could possibly be used as a genetic marker for wool trait selection in these sheep.

Key Grazing Behaviours of Beef Cattle Identify Specific Genotypes of the Glutamate Metabotropic Receptor 5 Gene (GRM5).

Genotype-phenotype associations between the bovine genome and grazing behaviours measured over time and across contexts have been reported in the past decade, with these suggesting the potential for genetic control over grazing personalities in beef cattle. From the large array of metrics used to describe grazing personality behaviours (GP-behaviours), it is still unclear which ones are linked to specific genes. Our prior observational study has reported associations and trends towards associations between genotypes of the glutamate metabotropic receptor 5 gene (GRM5) and four GP-behaviours, yet the unbalanced representation of GRM5 genotypes occurring in observational studies may have limited the ability to detect associations. Here, we applied a subsampling technique to create a genotypically-balanced dataset in a quasi-manipulative experiment with free ranging cows grazing in steep and rugged terrain of New Zealand's South Island. Using quadratic discriminant analysis, two combinations of eleven GP-behaviours (and a total of fifteen behaviours) were selected to build an exploration model and an elevation model, respectively. Both models achieved ∼ 86% accuracy in correctly discriminating cows' GRM5 genotypes with the training dataset, and the exploration model achieved 85% correct genotype prediction of cows from a testing dataset. Our study suggests a potential pleiotropic effect, with GRM5 controlling multiple grazing behaviours, and with implications for the grazing of steep and rugged grasslands. The study highlights the importance of grazing behavioural genetics in cattle and the potential use of GRM5 markers to select individuals with desired grazing personalities and built herds that collectively utilize steep and rugged rangelands sustainably.

Ovine KRTAP36-2: A New Keratin-Associated Protein Gene Related to Variation in Wool Yield.

Keratin-associated proteins (KAPs) are structural components of wool fibres. High-glycine/tyrosine (HGT)-KAPs are a subset of the KAP family, and their abundance in fibres varies. In this study, we report the discovery of an ovine HGT-KAP gene to which we assigned the name KRTAP36-2. Polymerase chain reaction and single-strand conformation polymorphism (PCR-SSCP) analyses revealed four variants of this gene in a screening population of 170 sheep from a variety of breeds. The DNA sequencing of the variants revealed four single-nucleotide polymorphisms (SNPs) and a dinucleotide deletion. Three of these SNPs were in the coding region, and one of these was non-synonymous and potentially led to the amino acid substitution p.Cys27Gly near the middle of the protein. The remaining SNP was located near the putative TATA box, and the di-nucleotide deletion was near the putative transcription initiation site. The effect of this variation in KRTAP36-2 was investigated in 274 Southdown × Merino lambs that were the progeny of five sires. Variation was only found to be associated with wool yield, that is, the proportion of the greasy fleece that remained as clean fleece upon scouring (expressed as a percentage). This may have some value in increasing wool production.

Variation in ovine KRTAP8-1 affects mean staple length and opacity of wool fiber.

In this study, keratin-associated proteins gene (KRTAP8-1) from five different sheep breeds and breed-crosses (n = 310) was genotyped using a Polymerase Chain Reaction-Single Strand confirmation Polymorphism (PCR-SSCP). Six unique genotypes were observed: AA, AC, AD, AE, DD and EE, with AA being the most common in the different breeds and crosses. Twelve wool characteristics: yield, mean staple length (MSL), bulk, mean fiber diameter (MFD), fiber diameter standard deviation (FDSD), coefficient of variation of fiber diameter (CVFD), medullation, standard deviation of medullation (MeSD), coefficient of variation of medullation (CVMed), opacity, standard deviation of opacity (OpSD), and coefficient of variation of opacity (CVOp) were measured on wool derived from the sheep. Variation in KRTAP8-1 was found to have strong association with MSL, OpSD and CVOp (p ≤ 0.027). The MSL of sheep of genotype AE was greater (p = 0.027) than for sheep of genotype AA. The OpSD of sheep of genotype AA was less (p = 0.017) than sheep with the AE genotype, and the CVOp of sheep with genotype AA was less (p = 0.018) than sheep with genotype AE. Further studies are required to confirm the role of variation in KRTAP8-1 in improving quality wool production.

Supplementing Diets with Agriophyllum squarrosum Reduced Blood Lipids, Enhanced Immunity and Anti-Inflammatory Capacities, and Mediated Lipid Metabolism in Tan Lambs.

Agriophyllum squarrosum (sand rice), a widespread desert plant, possesses anti-hyperglycemic and anti-inflammatory properties, and has been used in traditional Chinese medicine for many years. However, its effects on ruminants are unknown. To fill this gap, we examined the effects of A. squarrosum on the immune and anti-inflammatory responses of lambs. A total of 23, 6-month-old Tan ewe-lambs (27.6 ± 0.47 kg) were divided into four groups and offered a basic diet (C­control), or a diet that contained 10%, 20%, or 30% A. squarrosum, on a dry matter basis, for 128 days. Serum concentrations of total cholesterol were lower (p = 0.004) in the 30% supplemented lambs than controls, while concentrations of high-density lipoprotein cholesterol were lower (p = 0.006) in the 10% and 20%, but not in 30% supplemented lambs than controls. Serum-cortisol concentrations were lower (p = 0.012) in the 30% supplemented lambs and free fatty acid concentrations were higher in the 10% and 20% supplemented lambs than in control lambs (p < 0.001). Supplementation with A. squarrosum decreased (p < 0.05) the area of adipocytes in subcutaneous adipose tissue, but there was no difference between the 20% and 30% diets. Conversely, the area in visceral adipose tissue (VAT) increased (p < 0.05), especially for the 10% and 20% supplemented diets. Supplementation with A. squarrosum also enriched immune and anti-inflammatory related and lipid and glucose-metabolic pathways and associated differentially expressed gene expressions in adipose tissue. A total of 10 differential triacylglycerol, 34 differential phosphatidylcholines and seven differential phosphatidylethanolamines decreased in the diet with 30% supplementation, when compared to the other diets. Finally, adipocyte-differentiation genes, and immune and inflammatory response-related gene expression levels decreased in lamb adipocytes cultured with an aqueous A. squarrosum extract. In conclusion, supplementing lamb diets with A. squarrosum reduced blood lipids, enhanced immunity and anti-inflammatory capacities, and mediated lipid metabolism in adipose tissue and adipocytes of Tan lambs. A level of approximately 10% is recommended, but further research is required to determine the precise optimal level.

Variation in Ovine DGAT1 and Its Association with Carcass Muscle Traits in Southdown Sheep.

Diacylglycerol O-acyltransferase 1 (DGAT1) is a microsomal enzyme that plays a key role in the synthesis of triglycerides. Its gene (DGAT1) is regarded as a candidate gene for variation in milk and meat traits in cattle. The objective of this study was to use a PCR single-strand conformation polymorphism approach to explore sequence variation in two regions of ovine DGAT1 and to assess its effect on meat traits in New Zealand Southdown sheep. Three variant nucleotide sequences were identified in each region, with two single nucleotide polymorphisms (SNPs) and one nucleotide deletion being detected in intron 1 and two SNPs being found in exon 17. The effect of the exon 17 variation was not investigated due to one variant being predominant and the other two variants occurring at low frequencies. In intron 1, one variant (B1) was found to be associated with increase loin meat yield, suggesting that this may have value as a gene marker for improving meat traits.

The glutamate metabotropic receptor 5 (GRM5) gene is associated with beef cattle home range and movement tortuosity.

BACKGROUND: The grazing behaviour of herbivores and their grazing personalities might in part be determined genetically, but there are few studies in beef cattle illustrating this. In this study, we investigated for first time the genetic variation within a candidate 'grazing gene', the glutamate metabotropic receptor 5 gene (GRM5), and tested associations between variation in that gene and variation in grazing personality behaviours (GP-behaviours) displayed by free-ranging cows during winter grazing in the steep and rugged rangelands of New Zealand. Mature beef cows (n = 303, from 3 to 10 years of age) were tracked with global positioning system (GPS) and, with 5-minutes (min) relocation frequency, various GP-behaviours were calculated. These included horizontal and vertical distances travelled, mean elevation, elevation range, elevation gain, slope, home range and movement tortuosity, variously calculated using daily relocation trajectories with repeated measurements (i.e., 7 to 24 days (d)) and satellite-derived digital elevation models (DEM). The different GP-behaviours were fitted into mixed models to ascertain their associations with variant sequences and genotypes of GRM5. RESULTS: We discovered three GRM5 variants (A, B and C) and identified the six possible genotypes in the cattle studied. The mixed models revealed that A was significantly associated with elevation range, home range and movement tortuosity. Similarly, GRM5 genotypes were associated (P < 0.05) to home range and movement tortuosity, while trends suggesting association (P <  0.1) were also revealed for elevation range and horizontal distance travelled. Most GP-behaviour models were improved by correcting for cow age-class as a fixed factor. The analysis of GP-behaviours averaged per cow age-class suggests that grazing personality is fully established as beef cows reached 4 years of age. Home range and movement tortuosity were not only associated with GRM5 variation, but also negatively correlated with each other (r = - 0.27, P < 0.001). CONCLUSIONS: There seems to be a genetically determined trade-off between home range and movement tortuosity that may be useful in beef cattle breeding programmes aiming to improve the grazing distribution and utilisation of steep and rugged rangelands.

Sequence Variation in Caprine KRTAP6-2 Affects Cashmere Fiber Diameter.

Keratin-associated proteins (KAPs) are a structural component of cashmere fibers and in part determine fiber attributes. The gene encoding the high-glycine/tyrosine KAP6-2 (called KRTAP6-2) has been described in sheep, but it has not been identified goats. In this study, a 252-bp open reading frame with similarity to ovine KRTAP6-2 was found on goat chromosome 1, with its upstream and downstream flanking sequences are closely related with ovine KRTAP6-2 but are clearly distinct from other ovine KRTAP6-n sequences. Polymerase chain reaction amplification followed by single strand conformation polymorphism analysis of this region revealed five distinct banding patterns representing five different sequences (A to E) in 230 Longdong cashmere goats. Eleven diallelic single nucleotide polymorphisms (SNPs), a three-nucleotide sequence variation, and a 12-bp insertion/deletion were found among these five sequences, with most SNPs being either outside the coding region or synonymous. The presence of variant D was found to be associated with decreased mean fiber diameter (MFD; present: 13.26 ± 0.07 µm; absent: 13.55 ± 0.04 µm; p < 0.001), suggesting that variation in KRTAP6-2 may affect fiber diameter and have value as a molecular marker for improving the cashmere fiber diameter trait.

A Preliminary Investigation of Myostatin Gene (MSTN) Variation in Red Deer (Cervus elaphus) and Its Implications for Venison Production in New Zealand.

Myostatin (MSTN), also known as growth differentiation factor 8 (GDF-8), is a negative regulator of lean muscle tissue growth. Variation in the gene has been studied in many domesticated species, because of its potential to dramatically increase muscle mass. It has, however, not been investigated in red deer (Cervus elaphus). In this study, variation in MSTN intron 1 was investigated in 211 male New Zealand red deer, for which phenotypic measurements of M. Longissimus dorsi (eye muscle) (width, depth, and area, together with 12-month weight) were recorded. Two sequence variants (named A and B) differing by one nucleotide (c.373 + 224) were identified in the intron 1 region of the gene resulting in three genotypes (AA, AB, and BB; frequencies of 63.5%, 30.8%, and 5.7%, respectively), but no association between this variation and any of the quantitative measurements was detected. These results suggest that the deer MSTN is less variable than for other livestock species and that its activity may be controlled to maintain a size-growth equilibrium.

Identification of differentially expressed genes at different post-natal development stages of longissimus dorsi muscle in Tianzhu white yak.

The regulatory mechanisms controlling post-natal muscle development in the yak (Bos grunniens) are still largely unknown, yet the growth and development of muscle is a complex process that plays a crucial role in determining the yield and quality of an animal's meat. In this study, we performed a transcriptome analysis based on the RNA sequencing (RNA-Seq) of yak longissimus dorsi muscle tissue obtained from calves (6 months of age; 6 M), young adults (30 months of age; 30 M) and adult (54 months of age; 54 M) to identify which genes are differentially expressed and to investigate their temporal expression profiles. In total, 1788 differentially expressed genes (DEGs) (|log2FC| ≥ 1, P-adjusted < 0.05) were detected by pairwise comparisons between the different age groups. The expression levels of 10 of the DEGs were confirmed using reverse transcription-quantitative PCR (RT-qPCR), and the results were consistent with the transcriptome profile. A time-series expression profile analysis clustered the DEGs into four groups that could be divided into two classes (P < 0.05): class 1 profiles, which had up-regulated patterns of gene expression and class 2 profiles, which featured down-regulated patterns. Based on that cluster analysis, GO enrichment analysis revealed 1073, 127, and 184 terms as significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in the class 1 profiles, while 714, 66, and 206 terms were significantly enriched in BP, CC, and MF in the class 2 profiles. A KEGG pathway analysis revealed that DEGs from the class 1 profiles were enriched in 62 pathways, with the most enriched being the phosphoinositide 3-kinase (PI3K) - protein kinase B (Akt)-signaling pathway. The DEGs from the class 2 profiles were enriched in 16 pathways, of which forkhead box protein O (FoxO) - signaling was the most enriched. Taken together, these results provide insight into the mechanisms of skeletal muscle development, as well suggesting some potential genes of importance for yak meat production.

Variation in caprine KRTAP1-3 and its association with cashmere fibre diameter.

Keratin-associated proteins (KAPs) are components of cashmere fibres. The gene encoding the KAP1-3 protein (KRTAP1-3) has been described in goats, but little is known about sequence variation in this gene and if it affects cashmere fibre traits. In this study, we used a polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) technique to screen for nucleotide sequence variation in caprine KRTAP1-3 in 327 Longdong cashmere goats, then analysed association between the genetic variation that was revealed and some cashmere fibre traits. Six PCR-SSCP patterns representing six different variant sequences of KRTAP1-3 (named A to F) were revealed. Among these variant sequences, seven single nucleotide polymorphisms (SNPs) were detected, with two of them being non-synonymous. Goats with genotype AC had higher mean fibre diameter (MFD) than those with genotype AB (P < 0.001), while goats with genotype AB had higher MFD than those with AA (P < 0.001). The presence of C (P < 0.001) and B (P = 0.006) in a genotype was associated with increased MFD, and together this suggests that variation in caprine KRTAP1-3 affects the key fibre trait of MFD.

Ovine Toll-like Receptor 9 (TLR9) Gene Variation and Its Association with Flystrike Susceptibility.

Toll-like receptors (TLRs) are a family of proteins that play a role in innate immune responses by recognising pathogen-associated molecular patterns derived from various microbes. Of these receptors, TLR9 recognises bacterial and viral DNA containing unmethylated cytosine-phosphate-guanine (CpG) motifs, and variation in TLR9 has been associated with resistance to various infectious diseases. Flystrike is a problem affecting the sheep industry globally and the immune response of the sheep has been suggested as one factor that influences the response to the disease. In this study, variation in ovine TLR9 from 178 sheep with flystrike and 134 sheep without flystrike was investigated using a polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) approach. These sheep were collected from both commercial and stud farms throughout New Zealand and they were of 13 different breeds, cross-breds and composites. Four alleles of TLR9 were detected, including three previously identified alleles (*01, *02 and *03) and a new allele (*04). In total six single nucleotide polymorphisms (SNPs) were found. Of the three common alleles in the sheep studied, the presence of *03 was found to be associated with a reduced likelihood of flystrike being present (OR = 0.499, p = 0.024). This suggests that variation in ovine TLR9 may affect a sheep's response to flystrike, and thus the gene may have value as a genetic marker for improving resistance to the disease.

The Complexity of the Ovine and Caprine Keratin-Associated Protein Genes.

Sheep (Ovis aries) and goats (Capra hircus) have, for more than a millennia, been a source of fibres for human use, be it for use in clothing and furnishings, for insulation, for decorative and ceremonial purposes, or for combinations thereof. While use of these natural fibres has in some respects been superseded by the use of synthetic and plant-based fibres, increased accounting for the carbon and water footprint of these fibres is creating a re-emergence of interest in fibres derived from sheep and goats. The keratin-associated proteins (KAPs) are structural components of wool and hair fibres, where they form a matrix that cross-links with the keratin intermediate filaments (KIFs), the other main structural component of the fibres. Since the first report of a complete KAP protein sequence in the late 1960s, considerable effort has been made to identify the KAP proteins and their genes in mammals, and to ascertain how these genes and proteins control fibre growth and characteristics. This effort is ongoing, with more and more being understood about the structure and function of the genes. This review consolidates that knowledge and suggests future directions for research to further our understanding.

Sequence Variation in the Bovine Lipin-1 Gene (LPIN1) and Its Association with Milk Fat and Protein Contents in New Zealand Holstein-Friesian × Jersey (HF × J)-cross Dairy Cows.

Lipin-1 is known to play a regulatory role in tissues that function in lipid metabolism. In dairy cows, the lipin-1 gene (LPIN1) is highly expressed in the mammary gland, but its function in milk production is less understood. In this study, we used PCR-single strand conformation polymorphism analysis to investigate sequence variation in three regions of bovine LPIN1 in New Zealand Holstein-Friesian × Jersey (HF × J)-cross dairy cows, including part of the 5' non-coding region, the region containing the LPIN1ß-spliced exon, and the sixth coding exon that encodes the putative transcriptional activating domain of the protein. No variation was found in the LPIN1ß-spliced exon, but two sequence variants containing one single nucleotide polymorphism (SNP) were identified in the 5' non-coding region and four sequence variants containing four non-synonymous SNPs were identified in the sixth coding exon. Among the three common variants of the sixth coding exon, variant C was found to be associated with an increase in milk fat percentage (presence 4.96 ± 0.034% vs. absence 4.81 ± 0.050%; p = 0.006) and milk protein percentage (presence 4.09 ± 0.017% vs. absence 3.99 ± 0.025%; p = 0.001), but no associations (p > 0.01) were detected for milk yield. These results suggest that variation in LPIN1 affect the synthesis of fat and proteins in milk and has potential as a gene-marker to improve milk production traits.

MicroRNA-432 inhibits milk fat synthesis by targeting SCD and LPL in ovine mammary epithelial cells.

The microRNA (miR)-432 is differentially expressed in the mammary gland of two breeds of lactating sheep with different milk production traits, and between the non-lactating and peak-lactation periods, but there have been no reports describing the molecular mechanisms involved. In this study, the effect of miR-432 on the proliferation of ovine mammary epithelial cells (OMECs) and the target genes of miR-432 were investigated. The effects of miR-432 on the expression of the target genes and the content of triglycerides in the OMECs were also analyzed. Transfection with a miR-432 mimic was found using CCK8 and Edu assays, to inhibit the viability of OMECs and reduce the number of proliferated OMECs. In contrast, a miR-432 inhibitor had the opposite effect to the miR-432 mimic, and together these results suggest that miR-432 inhibits the proliferation of OMECs. A dual luciferase assay revealed that the genes for stearoyl-CoA desaturase (SCD) and lipoprotein lipase (LPL) are targeted by miR-432. The transfection of miR-432 mimic into OMECs resulted in decreases in the expression of SCD and LPL, and three other milk fat synthesis marker genes; FABP4, LPIN1 and ACACA. The mimic also decreased the content of triglycerides. The miR-432 inhibitor had the opposite effect to the mimic on the expression of these genes and the level of triglycerides. This is the first study to reveal the biological mechanisms by which miR-432 inhibits milk fat synthesis in sheep.

Variation in the Ovine Glycogen Synthase Kinase 3 Beta-Interaction Protein Gene (GSKIP) Affects Carcass and Growth Traits in Romney Sheep.

The glycogen synthase kinase 3 beta (GSK3ß)-interacting protein (encoded by the gene GSKIP) is a small A-kinase anchoring protein, which complexes with GSK3ßand protein kinase A (PKA) and acts synergistically with cAMP/PKA signaling to inhibit GSK3ß activity. The protein plays a role in regulating glycogen metabolism, protein synthesis, the cell cycle, and in regulating gene expression. In this study, PCR-single strand conformation polymorphism (PCR-SSCP) analyses were used to screen for variation in exon 1 and exon 2 of GSKIP in 840 New Zealand (NZ) Romney sheep. Two SSCP banding patterns representing two different nucleotide variants (A and B) were detected in an exon 1 region, whereas in an exon 2 region only one pattern was detected. Variants A and B of exon 1 had one non-synonymous nucleotide difference c.37A/G (p.Met13Val). The birthweight of sheep of genotype AA (5.9 ± 0.06 kg) was different (p = 0.023) to sheep of genotype AB (5.7 ± 0.06 kg) and BB (5.7 ± 0.06 kg). The hot carcass weight (HCW) of sheep of genotype AA (17.2 ± 0.22 kg) was different (p = 0.012) to sheep of genotype AB (17.6 ± 0.22 kg) and BB (18.0 ± 0.29 kg), and the fat depth at the 12th rib (V-GR) of sheep of genotype AA (7.7 ± 0.31 mm) was different (p = 0.016) to sheep of genotype AB (8.3 ± 0.30 mm) and BB (8.5 ± 0.39 mm). The results suggest that the c.37A/G substitution in ovine GSKIP may affect sheep growth and carcass traits.

Identification of sequence variation in the oocyte-derived bone morphogenetic protein 15 (BMP15) gene (BMP15) associated with litter size in New Zealand sheep (Ovis aries) breeds.

Advances in the study of reproductive traits indicate that functional variation in fertility genes may be useful for improving sheep fertility. The aim of this study was to search for variation in the bone morphogenetic protein 15 gene (BMP15) and ascertain any association with litter size in purebred Finnish Landrace sheep (n = 148), Finnish Landrace × Texel-cross sheep (n = 45), and composite sheep (of varying breed background; n = 58) from New Zealand (NZ). A 482 bp and 312 bp fragment of exon 1 and 2, respectively, of BMP15 were analysed using polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). The additive and dominance effect of BMP15 variation on litter size were estimated using animal and sire models. Two variants (A and B) were detected in exon 1; no sequence variation was detected in exon 2. Variant A had the nucleotide sequence CTT between positions c.31 and c.33, while variant B had a deletion (c.31_33del). The observed frequency for variant A in the Finnish Landrace sheep, Finnish Landrace × Texel-cross sheep and the composite sheep, was 0.77, 0.92, and 0.68, respectively while the frequency of variant B (c.31_33del) was 0.23, 0.08, and 0.32, respectively. An association between litter size and c.31_33del (P < 0.001) was observed in composite sheep. Analysis of more sheep will be required to confirm these results. Litter size did not differ significantly between sheep breeds regardless of the presence/absence of c.31_33del. Results suggested that c.31_33del might be a genetic marker for improving fecundity in some NZ sheep.

Variation in bovine leptin gene affects milk fatty acid composition in New Zealand Holstein Friesian × Jersey dairy cows.

Leptin is a protein hormone secreted from white adipose tissue. It regulates food/feed intake, body weight, immune function and reproduction. In our investigation, the polymerase chain reaction (PCR) amplification coupled with single-strand conformational polymorphism (SSCP) analysis was used to reveal variation in bovine leptin gene (LEP) in New Zealand (NZ) Holstein Friesian  × â€¯Jersey (HF  × â€¯J) dairy cows. Subsequent sequence analysis of a 430 bp amplicon spanning the entirety of exon 3 and part of the intron 2 region revealed three variant sequences ( A 3 , B 3 and C 3 ) containing a total of five nucleotide substitutions, all of which have been reported previously. Using general linear mixed-effect model analyses, the presence of variant A 3 (the most common variant) was associated with a decreased level of C15:1, C18:1 trans-11, C18:1 all trans, C18:2 trans-9, cis-12, C22:0 and C24:0 levels but increased levels of C12:1 and C13:0 iso ( p < 0.05 ). Variant B 3 was associated with reduced levels of C6:0, C8:0, C11:0, C13:0 and C20:0 but increased C17:0 iso and C24:0 levels ( p < 0.05 ). Variant C 3 was associated with decreased C17:0 iso levels but increased C20:0 ( p < 0.05 ) levels. In a genotype model, the A 3 B 3 genotype was associated with increased levels of C22:0 and C24:0 but decreased C8:0, C10:0, C11:0, C13:0, C15:0 and grouped medium-chain fatty acid (MCFA) levels ( p < 0.05 ). Genotype A 3 C 3 was found to be associated with decreased levels of C10:0, C11:0, C13:0 and grouped MCFA ( p < 0.05 ). This is the first report of findings of this kind in NZ HF  × â€¯J cows, and they suggest that variation in exon 3 of bovine leptin gene could be explored as a means of decreasing the concentration of saturated fatty acids in milk.

Nucleotide sequence variation of the major histocompatibility complex class II DQA1 gene in different cattle breeds from Nigeria and New Zealand.

The major histocompatibility complex (MHC) plays a role in immune response. Among other activities, the bovine MHC genes (BoLA) trigger immune responses, including the activation of antibody-producing B-cells. In this study, White Fulani (n = 24), Red Bororo (n = 5) and Holstein-White × Fulani-cross (n = 11) cattle from Nigeria, and New Zealand Holstein-Friesian × Jersey-cross (n = 40) cattle were used to investigate variability in exon 2 of BoLA-DQA1. Ten alleles were identified using a PCR-Single Strand Conformation Polymorphism (SSCP) approach and their nucleotide sequences confirmed by DNA sequencing. A total of 12.60 % of all nucleotide positions analysed were revealed to be variable and two novel BoLA-DQA1 alleles are reported here for the first time.